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PREFACE 


In  this  RAND  Memorandum  the  authors  present  further 
mathematical  results  In  the  problem  of  radiative  transfer 
In  a  one-dlmenslonal  medium.  This  subject  has  Important 
Implications  for  meteorology,  astrophysics,  and  the 
detection  of  nuc3ear  blasts. 

The  research  presented  here  was  sponsored  by  the 
Advanced  Research  Projects  Agency. 


SUMMARY 


In  the  present  Memorandvon,  by  means  of  the 
Invariant— Imbedding  technique,  the  Integral  equations 
for  the  reflection  and  transmission  coefficients  of 
radiation  In  a  one-dlmenslonal  medlmn  are  obtained, 
allowing  for  the  release  of  absorbed  energy  with  a 
random  time  delay.  Furthermore,  the  reflected  auid 
transmitted  Intensities  for  the  fluorescence  problem  in 
a  one— dimensional  medium  are  expressed  In  terms  of 
these  coefficients,  assuming  no  distribution  of  emitting 
sources  In  the  medium. 
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INVARIANT  IMBEDDING  AND  TIME-DEPENDENT  SCATTERING  OF 
LIGHT  IN  A  ONE-DIM12NSIONAL  MEDIUM 

1.  INTRODUCTION 

In  a  recent  series  of  papers,  based  on  the  principle 
of  Invariant  imbedding,  vdilch  stems  from  the  Invariance 
principles  of  Ambarzumian  and  Chandraselchar,  various 
kinds  of  time— dependent  neutron-transport  problems  in  a 
fixed  rod  of  finite  length  made  of  fissionable  material 
were  exactly  treated  by  Bellman  ai^d  Kalaba  [  l]  and  V/ing 
[ 12] .  It  was  shown  that  the  Laplace  transform  of  the 
Integral  equation  for  the  reflected  flux  derived  by  them 
provides  an  analytic  expression  of  the  solution  because 
of  the  convolution  form  of  the  integral  term  [  lo] . 
Furthermore,  the  technique  was  extended  [2]  to  the 
derivation  of  a  functional  equation  governing  the 
reflected  neutron  flux  from  a  rod  of  varying  length; 
the  einalytlcal  study  of  the  solution  of  the  partial 
differential  Integral  equation  was  made  with  the  aid  of 
an  iterative  procedure  [  111 . 

In  the  theory  of  radiative  transfer.  Miss  Busbridge 
[ 5I  used  the  Laplace  transform  for  reducing  the  non¬ 
stationary  transfer  equation  to  the  stationary  one.  The 
invariant— imbedding  technique  and  the  prlnclple-of— 
invariance  method  were  applied  by  Bellman,  Kaaaba,  and 
Ueno  [ 3I i  [  9^  to  time-dependent  diffuse  reflection 
problems  of  parallel  rays  by  an  Inhomogeneous  flat  layer. 
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On  the  other  hand,  time— dependent  scattering 
problems  of  light  In  a  one-dimensional  medivan  were 
discussed  by  several  Russian  astixjphysiclsts.  Intro¬ 
ducing  the  duration  of  temporal  capture  t^  and  the 
mean  free  time  t2  in  formulation  of  the  transfer 
equation,  Sobolev  [8]  extended  the  probabilistic  method 
to  some  transient  scattering  problems.  For  simplicity, 
he  derived  an  exact  solution  in  a  semi— infinite 
homogeneous  medium  when  t^  »  t^*  With  the  aid  of  the 
Laplace  transform  method,  Illnln  [  7^  obtained  an  explicit 
expression  of  the  quantvun  reflection  probability  and  of 
the  quantum  emergence  probability  at  any  level  from  the 
semi— infinite  homogeneous  medium,  allo\d.ng  for  both  time 
parameters,  t^^  and  tg*  Recently,  Kaplan  [6]  also 
extended  Sobolev's  procedure  to  the  three-dimensional 
nonsteady  case  eind  discussed  the  functional  equation  for 
the  reflection  coefficients  in  a  semi— infinite  homogeneous 
medium  vAien  t^^  »  tg* 

In  a  series  of  papers,  making  use  of  the  probabilistic 
method,  Kaplan  and  others  [63  considered  the  scattering 
of  light  in  a  one— dimensional  semi- infinite  homogeneous 
medltun  with  a  moving  boundary.  Blberman  and  Veklenko  [4] 
also  extended  their  probabilistic  approach  to  the 
derivation  of  relationships  for  the  generalized 
reciprocity  principle  for  noncoherent  scattering  in 
radiative  transfer.  Lists  of  many  other  papers  attempting 
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to  solve  the  time— dependent  transfer  process  will  be 
found  in  the  references  of  our  precedinc  papers  L >  19^* 
In  the  present  Memorandum,  the  invariant— imbedding 
technique  is  used  to  derive  the  integral  equations  for 
the  coefficients  of  reflection  and  transmission  of 
radiation  by  a  one— dimensional  homogeneous  medium  of 
finite  optical  thiclaiess,  allov/ing  for  tvro  characteristic 
time— scales,  t^  and  t^.  Furthermore,  the  reflected 
and  transmitted  intensities  from  this  medium,  illuminated 
by  a  imlt  step— function  time-dependent  pencil  of 
radiation  from  outside,  are  expressed  in  terms  of  these 
coefficients,  assuming  no  distribution  of  emitting 
sources  in  the  medium. 

2.  THE  EQUATION  Ql<'  TRANSFER 

Consider  a  one— dimensional  homogeneous  medium  of 
optical  thickness  >  0,  illuminated  by  radiation  of 
time-dependent  specific  intensity  1^  Incident  on  the 
right-hand  boundary  t  =  (see  Fig.  1).  Scattering 
of  light  in  either  direction  is  assumed  to  be  equally 
probable. 

,  ,  Ip(T,t)  I,(T,t)  Il(T,,t) 

ipCOit)  J ^  ^  ^  , 

^ -  I - 1 - 1  - - 

0  T  l2(Tl>t) 

Fig.  1 

A  Time-dependent  Transport  Process 
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Lc-t  I^{T,t)  and  denote  respectively  the 

specific  intensities  of  radiation  at  the  level  t  at 
tine  t,  dli’ected  tov;ard  the  boundaries  t  =  and 
T  =  0. 

Tlie  equation  of  transfer  appropriate  to  the  present 
case  takes  the  form 

dl  ,  c)I 

(2*^)  ZPT"  ^  jF"  +  ^1  " 

*^'1  I 

(2.2)  +  B(T,t), 


vriiere  c  is  the  speed  of  llcht,  i  is  the  attenuation 
coefficient,  and  3(T,t)  is  the  source  function. 

V/e  assume  that  the  duration  of  temporal  capture, 
i^iich  corresponds  to  the  mean  molecular  interaction  time 
in  the  kinetic  theory-  of  dilute  cases,  is  equal  to  t^, 
and  the  probability  of  emission  during  the  Interval  of 
time  (t,t  +  dt)  is  clven  by  exp(—  t/tj^)dt/t^. 

Assiimlng  the  above  emission  probability  in  the 
successive  scattering  processes,  consisting  of  absorption 
and  emission,  we  see  that  the  source  function  B(T,t) 
can  be  witten  in  the  fom 


(2.3)  B(T,t)=|J  [l^(T,t')  +  l2(T,t')]e  ^ 

-00  ^ 


vrfiere  a  is  the  albedo  for  single  scattering,  l.e.,  the 
probability  of  quantum  survival.  Under  some  limited 


conditions  Sobolev  [8]  discussed  the  solution  of 

(2.1)  -  (2.2). 

If  we  v/rite 

(2.4)  ci  =  1/t^j 

nhei'e  t2  is  the  mean  free  time,  the  equations  of 
transfer  (2.1)  —  (2.2)  become 

dl 

(^•5)  ^2  TPT  =  -  +  B(T,t), 

I  "J I 

(2.6)  —  -y:—  +  t2  -pp^  =  —  I2  +  3(T,t), 

tonether  .with  the  boundary  and  initial  conditions 

(2.7)  1^(0, t)  =  0,  12(^1, t)  =  1^(4), 

(2.8)  I^(t,0)  =  0,  I2(T,0)  =  0  (0<T^T^), 

where  called  respectively 

the  reflected  and  transmitted  Intensities,  and  l2(T^,t) 
is  the  intensity  incident  on  the  boundary  t  =  at 
time  t. 


3.  THE  REFLECTED  INTMSITY 

Let  R(T^,t)  denote  the  coefficient  of  reflection. 
Consider  that  at  time  t  =  tQ  (0  tQ  ^  t)  the  incident 
radiation  falls  on  the  boundary  t  *  t^.  Because  the 
optical  properties  of  the  medium  are  constant  with 
respect  to  time,  for  convenience  we  can  put  tQ  equal 
to  zero  without  loss  of  generality.  Then,  by  definition. 


vrtiere  l2(T^,t)  Is  given  by  (2.7). 

We  shall  seek  an  integral  equation  for  the 
reflection  coefficient  R(T2,t),  making  use  of  the 
invariant— imbedding  technique. 

Imbedding  the  one-dimensional  medixan  of  the  optical 
thlclmess  in  position  and  time,  we  have 

(3.2)  I^(t^  +  A,t  +  tgA) 

=  Ij_(T^,t)  +  A.{-  +  B(T^,t))  +  o(A), 

where  o(A)  is  of  the  order  of  magnitude  of  the 

p 

infinitesimal  A  ,  From  (2.6)  we  get 
(3*3)  ”  ^2^^ 


-  l2(Tj^,t)A  +  B(T^,t)A  +  o(A), 
By  initial  condition  (2.8),  (3.3)  becomes 
(3A)  l2(T^it)  =  lQ(t  -  t2A)  -  I^(t)A 


+  f 


-(t-t')/t, 

Wt.)e 

-00 

r‘  ,  , at. 


Making  use  of  (3.1)  and  (3.^)#  we  obtain 


r  ^  '' 

(3.5)  "  J 

—00 

“  I  “  t')lo(t‘  -  t2A)dt' 

-00 

-  R(T;^,t  -  t')lo(t')dt' 

-00 

+  I  A  J  R(T^,t  -  t' )dt' 

-00  ^ 

+  I  A  J'  ^  R(T^it  -  t' )dt» 

—00 


—00 


On  the  other  hand,  we  have 


(3*6)  ^2^^ 
p  t+tpA 

“j  R(*f3^  +  A,t  +  tgh  -  t')lo(t')dt*. 


-00 


With  the  aid  of  (3.5)  and  (3.6),  (3.2)  becomes 


-8- 


p  t+toA 

(3.7)  J  ^  RCTj^  +  A,t  +  tgA  -  t')lo(t')dt» 
—00 
t 


=  f  R(T^,t  -  t«)lQ(t'  -  t2A)dt' 

—00 

r  ^ 

-2a  J  R(T^,t  -  f  )lo(t')dt' 

-00 

rt  ft'  -(t'-t")/t, 

+  |aJ  R(T^,t  -  t')dt'  J  I^(T^,t")e  ^ 

-00  -00  ^ 
ft  ft'  -(t'-t")/t,  ^+.1- 

+  |aJ  R(T3,t  -  t')dt'  j 


-00 
t 


+  fAj  Ii(T^,t')e 

-00 


-(t-t')/t^  dt' 


-(t-t')/tj^  dt' 


-00 


Hie  Dirac  delta— function  time-dependent  case.  In 
v/hat  follov^s,  \re  shall  treat  the  case  of  the  Dirac  delta 
time— dependent  function,  lQ(t)  “  where  F  is  a 

constant  and  6  is  the  Dirac  delta— function.  The 
substitution  of  lQ(t)  =  F6(t)  into  (3.7)  provides 


(3.8) 


R(t^  +  A,t  +  tgA)  =  R(T^,t  -  tgA)  -  2AR(T^,t) 

a  .  r^'  ~(t'-t")/t, 

+  |Aj  R(T^,t  -  t')dt'  J  R(T^,t")e 

-00  -00  ^ 

+  I  aJ  ^  R(T^,t  -  tOe"*  ^ 

.f.r 
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Hence,  letting  A  -♦  0,  we  get 


(3.9) 


Sr  .  0+.  ^R  A  OR 
BtT  +  2t2  3^  +  2R 


a 


— t/t^ 


TEr~+J  R{‘fi,t  -  t  )e 

•*■  -00  ■*■ 


-f/t. 


,  pt  pt'  -(t'-t")/t, 

+  |j  dt'  J  R(Tj^,t  -  t*  )R(T3^»t")e 


-00  -00 
Hie  conditions  Inyiosed  on  R  are 


(3.10)  R(T^,t)  »  0  for  0  >  t;  R(0,t)  »  0  for  t  >  0. 


Equation  (3.9)  Is  the  requisite  Integral  equation 
governing  the  reflection  coefficient  R(Tj^,t). 


The  unit  step-function  time-dependent  case.  Consider 
a  fluorescence  problem  for  which  the  diffusely  reflected 
light  decreases  for  a  long  time  after  the  sudden 
switching— off  of  the  external  radiation  Incident  on  the 
boundary,  assvonlng  no  emitting  sources  In  the  medium.  In 
the  time  Intervea  (-  oo,  0),  let  the  radiation  of 
constant  Intensity  F  fall  on  the  boundary  t  =  t^,  and 
at  time  t  »  0  let  the  Incident  radiation  be  suddenly 
quenched. 

Writing 

(3.11)  13(^1, t)  -  PH*(t), 


where 


-2  y 


(3.12) 


*  ,  V 

H  (t) 


1,  t  <  0, 
0,  t  >  0, 


we  find  that  the  requisite  Intensity 

reflected  by  the  end  t  =  at  time  t,  is  given  by 

(3.13)  Ii(T;^it)  “  R(T^,t-t«)H*(t')dt' 

-00 
n  00 

=  F  J  R(Tj^,u)du, 
t 


where  R2('r]^jt)  is  governed  by  (3.9). 

In  the  next  place,  consider  a  Heaviside  unit-function 
time~dependent  case,  in  tdiich  from  the  Instant  t  «  0 
the  radiation  of  constant  intensity  F  falls  continuously 
on  the  boundary  t  ■  from  the  outside.  We  ask  for 
the  gradual  increase  of  the  intensity  reflected  from  the 
medium  at  time  t  >  0.  Writing 


(3.14) 
where 

(3.15) 


12(^1, t)  - 


H(t) 


FH(t), 

t  >  0, 

t  <  0, 


we  obtain 


-  R(Tj,t-t')H(t')dt» 

-00 

■  F  J'  R(Tj^,u)du, 

0 


(3.16)  Ii(T^,t) 
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whlch  provides  the  reqxilred  reflected  Intensity  in  this 
case. 

4.  THE  TRANSMITTED  INTENSITY 

Let  T(T^,t)  denote  the  coefficient  of  transmission. 

Then  we  have 

(4.1)  l2(0,t)  T(T^,t-t')l2('fi>t')dt'. 

—00 

We  inqiiire  into  an  integrail  equation  for  T(T^,t). 

In  a  manner  similar  to  that  used  in  the  preceding  section, 
we  have 

(4.2)  I2(0,t+t2d)  =  l2(0,t)  +  A{-  l2(0,t)  +  B(0,t))  +  o(d). 
From  (4.1),  we  obtain 

n  t+tpd 

(4.3)  I2(0,t+t2d)  =  j  T(  ^+d,t+t2A-t')lo(4')cit'. 

—00 

On  the  other  hand,  using  (3.4),  we  see  that  the  trans¬ 
mitted  intensity  l2(0,t)  is  given  by 
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(4.4)  l2(0,t)  T(T^,t-t»)l2(T^,t')dt' 

-00 

T(Ti,t>-t')lo(t'-t2A)dt' 

-00 

-  T{T^,t-t')lo(t')dt' 


+  f  A  y'  T(T^,t-t')dt 
-00 

ft*  -(t'-t")/t, 

-00 

+  ^  aJ^  T(Tj^,t-t'  )dt' 

-00 


Then,  recalling  (2.8),  (4.3),  and  (4.4),  we  find  that 
(4.2)  becomes 


(4.5)  / 


t+t2A 


T( Tj4A, t+t2A-t • ) Iq( t » ) dt » 

-00 

f  T(T^,t-t«)lQ(t'-t2A)dt'  -  T(T^,t-t')lQ(t')dt' 

-^0  ~^x> 

+  T(Ti,t-t')df  Ii(T^,t")e''^*'"^  ^/^l  ^ 

-00  -00  ^ 

+  I  a  J  ^  T(Tj,t-t»)dt«  lQ(t")e~^^'~^  ^ 

-00  -00  ^ 

+  a|-  l2(0.t)  +f/‘  l2(0,f  )•”'*■*' 3|l|  +  0(4). 

I  1 
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The  Dirac  delta-function  time-dependent  case. 
Inserting  -  P6(t)  Into  (4.5)  and  letting 

A  -»  0,  we  have 


-00  ^ 
+  |J  dt'  J  T{T^,1>-t')R(T^>t")e 


—00  —00 


along  with  the  boundary  and  initial  conditions 


(4.7)  T(T^,t)  =0  (t^  >  0,  0  >  t  or  t  < 

T(0,t)  =  Fa(t)  (t  >0). 


The  vinit  step-fvinction  time-dependent  case.  To 
begin  with,  we  consider  the  same  quenching  fluorescence 
problem  as  that  treated  in  the  preceding  section. 

Under  the  incident  intensity  given  by 

(3.11),  the  requisite  Intensity  transmitted  from  the 
boundary  t  =  0  at  time  t  is  provided  by 

ft  * 

(4.8)  1^(0, t)  =  F  j  T(T^,t-t')H  (t')dt' 

-00 
r  00 

=  P  j  T(T^,u)du, 
t 

when  T(T^,t)  satisfies  (4.6). 

Next,  we  seek  the  transmitted  intensity  when  the 
boundary  t  =  is  Illuminated  by  a  pencil  of  external 
radiation  of  intensity  l2(Tj^jt),  given  by  (3.13). 
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Then  we  have 

(4.9)  13(0,4)  =  P/^  P(Tj^,1>-t')H(t')dt' 

-00 

»  F  y'  T{T2^,u)du. 

0 

5.  DISCUSSION 

In  later  papers,  the  present  approach  will  be 
applied  to  problems  of  light  scattering  in  a  finite 
one-dimensional  medium  with  a  moving  boundary. 
Purthemore,  the  fluorescence  problem  in  vSiich  the 
distribution  of  emitting  sources  in  the  one-dimensional 
medium  is  a  function  of  t  and  the  external  light 
source  is  distinguished  at  time  zero  will  be  considered. 
An  analytical  and  computational  study  of  the  solutions 
of  the  integral  equations  for  the  R-  and  T-functlons 
will  also  be  made. 
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